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Quantized Gauged Massless Rarita-Schwinger Fields 

Stephen L. Adle 

Institute for Advanced Study, Einstein Drive, Princeton, NJ OSSfO, USA. 

We study quantization of a minimally gauged massless Rarita-Schwinger field, by both 
Dirac bracket and functional integral methods. The Dirac bracket approach in covariant 
radiation gauge leads to an anticommutator that has a non-singular limit as gauge fields 
approach zero, is manifestly positive semidefinite, and is Lorentz invariant. The constraints 
also have the form needed to apply the Faddeev-Popov method for deriving a functional 
integral, using the same constrained Hamiltonian and inverse constraint matrix that appear 
in the Dirac bracket approach. 


I. INTRODUCTION 


In this paper we continue the study of gauging a massless Rarita-Schwinger field begun in 
the preceding paper [l|, referred to henceforth as (I) , which dealt principally with the classical 
case (with a small excursion into first quantization). Here we turn to a detailed examination of 
quantization of a gauged massless Rarita-Schwinger field. Our main aim is to show that a consistent 
quantization is possible in gauge covariant radiation gauge, avoiding the problem of non-positivity 
of the canonical anticommutator first noted by Johnson and Sudarshan ^ and later rederived by 
Velo and Zwanziger Other objections to gauging a massless Rarita-Schwinger field - the issue 
of superluminal signaling, and various “on-shell no-go” theorems - have already been taken up in 
(I). In referring to a formula numbered “Eq. (^)” in the preceding paper we shall use the notation 
“Eq. {l-ffy , while non-hyphenated equation numbers refer to equations from this paper. 

In Sec. 2 we give the Hamiltonian form of the equations of motion and constraints, and introduce 
the Dirac bracket. This can be done without imposing a gauge fixing condition; in articular, 
we do not use the condition Tq = 0 that was imposed in an initial arXiv posting ^ of this 
paper. When a gauge hxing condition is omitted, the equation of motion for T computed from the 
Dirac bracket agrees with the equation of motion of Eq. (T29), in the form obtained when Tq is 
eliminated by using the secondary constraint w = 0; this demonstrates that the Dirac formalism is 
working correctly in the Rarita-Schwinger equation context. However, in the absence of a gauge 
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fixing constraint, the Dirac bracket anticommutator of 'h with agrees with the anticommutator 

n n 

calculated in [3|] and [^, which is singular in the limit of vanishing gauge fields and is not positive 
semidefinite. 

In Sec. 3 we study the Dirac bracket in its classical and quantum forms with imposition of 
a covariant radiation gauge constraint. We show that now the quantum Dirac bracket has the 
requisite positivity properties to be an anticommutator; related details are given in Appendix A. 
In Sec. 4 we give an alternative approach to proving positivity of the anticommutator in covariant 
radiation gauge, based on writing a Lagrangian for the equation of motion for 4' in which 4^0 has 
already been eliminated by use of the secondary constraint. In Sec. 5 we discuss Lorentz covariance 
of covariant radiation gauge and show Lorentz invariance of the Dirac bracket. In Sec. 6 we turn to 
path integral quantization in covariant radiation gauge, leading to a formalism closely resembling 
the Dirac bracket approach. A brief concluding discussion is given in Sec. 7. 

Our conclusion from this paper and the preceding one is that one can consistently gauge a 
massless Rarita-Schwinger field, at both the classical and quantum levels. This opens the possibility 
of using gauged Rarita-Schwinger fields as part of the anomaly cancelation mechanism in grand 
unified models, with anomalies of the spin ^ fields canceling against the spin | anomaly. 


II. HAMILTONIAN FORM OF THE EQUATIONS AND THE DIRAC BRACKET 

The standard route to canonical quantization is to transform the Lagrangian equations to 
Hamiltonian form, and to take the constraints into account by replacing the classical brackets by 
Dirac brackets. In carrying this out, we will simplify the formulas by making the gauge choice 
Aq = 0 for the non-Abelian gauge fields. This gauge choice is always attainable, and leaves a 
residual non-Abelian gauge invariance with time-independent gauge parameter. The Hamiltonian 
will then be covariant with respect to this restricted gauge transformation. For the moment, in 
discussing the canonical Hamiltonian and bracket formalism, we will allow A to be time dependent, 
so that E 0. But when we turn to the Dirac bracket construction corresponding to a constrained 
Hamiltonian, which is simplest in the case of time-independent constraints, we will assume a time- 
independent A, corresponding in Aq = 0 gauge to S = 0. (If we carry along the Aq term in 
the formulas then time-independent fields would not require E = 0. So this specialization can be 
avoided at the price of somewhat lengthier equations.) 

From the action S'(4'^) = f dtL {^of Eq. (1-23) and the canonical momentum P = x <?, 
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we find the canonical Hamiltonian to be 


H = / d?xdn^ -P-L 


- [ • Z5 X ^ • a X D^o + • Z5 x 


= -- / cPx[-%a ■ D X ^ + {iP - P X a) ■ {a X D^q + D x ^)] , 


where in the final line we have used the inversion formula = iP — P x a. 

We can now compute the classical brackets of various quantities with H. From 
^ 1 

— =[^,H]c = -[i{a X + Dx^)-ax{dx D^q + D x' if)] 
dt 2 

=Difo + ^[-1? X (z5 X ^) + X , 


( 1 ) 


( 2 ) 


we obtain the \I' equation of motion in the form given in Eq. (1-29). Similarly, from the bracket of 
P with H we find the equation of motion for Turning to brackets of the constraints with H, 
starting with P^f , we find 


dP, 


dt' 2 ^ ’ 


(3) 


and so P^i = 0 for all times implies that x = 0. For the total time derivative of X; we have 

^ = ^ + lx.H]c = s->^g^-i + lx,H]c = -igu, , 


(4) 


and so X = 0 for all times implies that oj defined in Eq. (T28) vanishes. Since oj contains a term 
proportional to Tq, to continue this process by calculating the time derivative of oj, we must obtain 
difo/dt from a bracket of Tq with H (and similarly for d^\^/dt). This requires adding to H a term 


AH = - d^x 




(5) 


Requiring AH to be self-adjoint then imposes the requirement 


Pi. = -P, 


'I'o 




( 6 ) 


which was noted following Eq. (1-56). As noted in (I), the chain of successive brackets with H 
starting from P^ t and continuing to X;*^; ••• leads only to constraints involving T and Tq but never 
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their adjoints. The doubling of the set of constraints, which turns the first class constraints into 
second class ones, comes from requiring that the adjoint of each fermionic constraint also be a 
constraint, not from taking successive brackets with H. 

We are now ready to implement the Dirac bracket procedure. The basic idea is to change the 
canonical bracket [F, G]c to a modified bracket [F, G\d-, which projects F and G onto the subspace 
obeying the constraints, so that the constraints are built into the brackets, or after quantization, 
into the canonical anticommutators. The constraints can then be “strongly” implemented in the 
Hamiltonian by setting terms proportional to the constraints to zero. After integration by parts 
the second line of Eq. ([T]) takes the form 

H = -]^j Sx[-^lx-X^^o + ^^-Dx^] , (7) 

so setting the constraints xK X respectively to zero in Eq. ([7]), we see that the constrained Hamil¬ 
tonian is just 

H = -^ j • F X ^ 

~ ~ ^ J d^x{iP — F X (t) • F X T 

( 8 ) 

which coincides with the energy integral computed in Eq. (T35) from the stress-energy tensor. 

We proceed now to calculate the Dirac bracket for the case when F = F('k) and G = G(T, '&!); 
the case when F = F('k^) can then be obtained by taking the adjoint, and the case when F = 
F('k,T^) can be obtained by combining the extra bracket terms from both calculations. When 
F has no dependence on it has vanishing brackets with the constraints (l)a of Eq. (1-55) and 
nonvanishing brackets with the constraints Xa of Eq. (1-56). The Dirac bracket then has the form 
(see Eqs. (I-A20) and (I-A21) for why M~^ appears) 

lF.Ob = |F,G|c-EE[^. Xa]cM-,^[^b,G] , (9) 

a b 

where Mab{x,y) = [(t’aix), Xb{y)]c is the matrix defined in Eqs. (1-58) and (1-59). We recall that 
this matrix has the form 

/0 -1 0 0 \ 


M = 


\ U S T 
Q V A B 
yO W C Fy 


( 10 ) 
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where in the SU{n) gauge field case, each entry in M is a 2n x 2n matrix. Using the block inversion 
method given in Eqs. (I-A18) and (I-A19), we find that M~^ is given by 


^ S 1 

- 1 0 
ri2 + gw 0 
\nv + iw 0 


-{ST + TH) -{Sg + TI)^ 
0 0 

j' g 

n z I 


(11) 


where 


Y. = u-s{Tv ^gyv)-T{nv ^zw) , 


( 12 ) 


and where T^ g, Ti, Z are the elements of the block inversion of the matrix N of Eq. (1-60), 




(13) 


Substituting these into Eq. ([9]) we find for the Dirac bracket a lengthy expression, which sim¬ 
plifies considerably after noting that [E(^),xi]c = [E(^), —P,j,p]c = 0 and [(^i, = 

, G'(^, ^l)]c = 0, leaving the relatively simple formula 


-[F{^),X3]c{z^ ^^)]c + g [04, G(^, 

-[F{^),XA]c{n [03, , ^t)]c + Z [04, G(^, ^t)]c) . 

(14) 


We note that only the matrix N enters, in this case through its inverse, rather than the full matrix 
of constraint brackets M. The final step is to evaluate the inverse block matrix elements P, g, H, Z 
from the expressions for A, B, C, V, again by using the block inversion formulas of Eqs. (I-A18) 
and (I-A19). Let us define the Green’s function Z)~^{x — y) by 

(i{L^f+ a ■L^xL^)V-^{x-y) = 5‘^{x-y) , (15) 

and a second Green’s function Z{x — y) by 


Z(f-y) =A-BV-^C 

= - 2iga ■ B6^ {x - y) - ■ LsV~^{x -y)L^-^^ 


(16) 
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where in covariant radiation gauge L = D. Then the needed inverse block matrices are 

T=Z-^ , 

Q = -Z-^BV-^ , 

n = -v-^cz-^ , 

X=V-^ ^V-^CZ-^BV-^ . 

(17) 

We wish now to apply the Dirac bracket formula to the cases (i) T(T) = T and G(T, 
and (ii) -F(T) = 'k and G('k,'k^) = 77, with H the constrained Hamiltonian of Eq. ([8]). The 
following canonical brackets are needed for this: 

[^(®),X3(2/)]c =2Di-<5^(T-y) , 

[^(^),X4(y)]c - o' X Ls)5^{x-y) , 

[(/> 3 (f),^^(y)]c =2Dg6^{x-y) = -26^{x-y)^^ , 

[4>A{x),'^\y)]c = - {iLs- Ls; X d)d^{x-ip) = 6^{x - y){i^^g x a) , 

[(t> 2 .{x), H]c =igB{x) ■ T(f) , 

[(/>4(x),i7]c X <t) X 75^-T(f) . 

(18) 


Additionally, for case (i) we need the canonical bracket 

[Tj(x),T](y)]c =[Tj(x),zPj(y) - ejkiPk{y)(yi\c 

= - i{Sij + iejiiai)6^{x - y) = -iajai6^{x - y) = -2i(^6ij - ^aiaj'jd^ix - y) , 

(19) 


and for case (ii) we need the canonical bracket 

1 


[^i{x),H]c = -^{iDx X ^(x) - (T X (ps X ^(x))^^ 


( 20 ) 


Up to this point, we have not specialized L so as to make it easy to ascertain what the formulas 


become when gauge fixing is omitted (as in 


Q] and [ 4 I). 


When L = 0, the matrix N degenerates 


to its upper left element A. This is reflected in the fact that Z of Eq. (1160 simplifies to 


Z{x — y) = A = —2igd ■ BS^{x — y) 


( 21 ) 
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which is a local function of x and so is algebraically invertible. The Dirac bracket of T(x) with the 
constrained Hamiltonian now simplifies to 


d^{x) 

dt 


=[^{x),H]d = -[iDg X T(x) - a X [Dg x T(x))] - J d^y^2Dg Z ^{x 
= 7 ;[iDg X T(f) - a X [Dg x T(f))] + Dg ^ \^^ B{x) ■ T(x) 


a ■ B{x) 


y)igB{y) ■ T(y) 


} 


( 22 ) 

The second line of this equation is just the T equation of motion in the form of Eq. (1-29) 
(when Hq = 0), with Tq eliminated by using the secondary constraint, which when E = 0 reads 
d ■ B^q = H • T. This shows that the Dirac bracket formalism correctly incorporates the Tq term 
of Eq. (T29). The reason a local result is obtained from this calculation is that in the absence of 
gauge fixing, the Dirac bracket only projects into the subspace that preserves the primary constraint 
X = 0, and since the equation of motion of Eq. (T29) preserves this constraint, it already resides 
in the subspace projected into by the L = 0 Dirac bracket. 

When L = 0, for the Dirac bracket of Ti(x) with Tj(y) we find 


[^iix),^]{y)]D =[^iix),^]{y)]c - J d^wd^z[^i{x),X3{w)]cZ ^(zu - T)[(()3(T), 


c 


= -2i 
= — 2i{x 


(^ij - ^(Xiaj'j6^{x-y) - 


T -.H^' 

2 V gd-B 


gd ■ B{x) 

1^ , 


(23) 


where in the final line we have written iDgi = Hj to relate to the abstract operator notation 
of Velo and Zwanziger Multiplying the final line by i to convert the Dirac bracket to an 
anticommutator, and by a factor 1/2 reflecting our different field normalization, Eq. (j23P becomes 
the expression for the anticommutator given in the zero mass limit of Eq. (4.12) of [^. Using 
identities in Appendix A of (I), one can verify (as in Appendix C of P]) that 

1 


{d X Dg)i [6ij - ^aiaj]S^{x - y) - Dgi ^ 


yj 


- RIT - ’ (24) 

ga ■ B[x) 

that is, the constraint x is explicitly projected to zero. However, as noted in the Introduction to 
(I), the anticommutator of Eq. ()23p becomes singular as B ^ 0, rather than limiting to the free 
Rarita-Scwhinger anticommutator. This problem is a direct consequence of omitting a gauge-fixing 
constraint, by taking L = 0 in calculating the matrix N. 
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Now setting L = D for covariant radiation gauge, we find for the Dirac bracket of 'l'i(ai) with 
the constrained Hamiltonian, 

=[^{x),H]d = liiD^ X ^(x) - a X {D^ x ^(x))] 

-/ d^y[‘^Ds; T{x-y)igB{y)--^{y) 

+G{x - y)^{iD^ - X a) X D^- ^{y) 

+ (zZ5^ - a xL^) nix - y)igB{y) ■ ^ (y) 

+I{x-y)^{iD^-D^x a) X D^-^{y) ^ . (25) 

The first line of this equation gives the second term of the unconstrained equation of motion in 
the form of Eq. (T29), while the remaining terms replace the first term of Eq. (T29) to guarantee 
that 

#3 ^ d{a xD ■ T) ^ ^ - n 

dt dt dt dt ’ 

dc^A _,D • ^ ^ _ n 

dt dt ' dt 

(26) 


where we have used the fact that we are assuming that D is time independent. That is, the Dirac 
bracket simultaneously projects the equation of motion into the subspace where both x = 0 and 
D ■ ^ = 0. The restriction to D time independent can be avoided by treating the gauge helds 
as dynamical variables, taking into account their own constraint structure, and noting that the 
radiation gauge fixing constraint V -P^ = 0, with the canonical momentum conjugate to A, has 
nonvanishing fermionic brackets with all Rarita-Schwinger constraints involving D = V + gA. This 
requires an extension of the Dirac bracket construction to take the new, Grassmann-odd, brackets 
into account, and the extended Dirac bracket structure will then obey Eq. (I26p without requiring 
the assumption of a time independent A and D. 

With L = D, we find for the Dirac bracket of 'l'i(x) with Tj(y), 

[Tj(x),Tj(y)]D = - - ^cJiCJj)(5^(x - y) 

T4D,j? iP (x y')^5yj ^BxiG{,x y) y y ^ *^)j 

-{-‘2{iBx (T X Bx')i'hL(^x y^^B y j (iD^ ex B^iP{x y^{iB g ^ 5 

(27) 


which gives the generalization of Eq. (I23p to the case when a covariant gauge fixing constraint is 
imposed. This equation will be further analyzed in the next section. 
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III. QUANTIZATION OF THE ANTICOMMUTATOR DERIVED FROM THE DIRAC 
BRACKET AND POSITIVITY IN COVARIANT RADIATION GAUGE 


Given the Dirac bracket, the next step is to quantize, by multiplying all Dirac brackets by i and 
then reinterpreting them as anticommutators or commutators of operators. In the case considered 
here, this can be done in a constructive way, as follows. First let us replace the set of 2n component 
column vector constraints (/iq and 2n component row vector constraints Xa by the set of 4n scalars 
given by their individual matrix elements. Moreover, since the Xa are the adjoints of the (pa, we can 
take linear combinations to make all of these scalars self-adjoint. Labeling the set of self-adjoint 
scalar constraints by <ba, the Dirac bracket construction for the bracket of F with G reads 


[F, G]d =[F, G]c - E ^a]cT-,^[^b, G]c 

a b 


Tab=[^a,^b]c , 


(28) 


with the matrix T real. 

We now observe that since the are all linear in the scalar components of ^ and if we 
make the replacement z[ , ]c { , }c, with { , } the anticommutatior, and replace all Grassmann 
variables ^ and with operator variables having the standard canonical anticommutators, then 
since there is no other operator structure the same real matrix Tab will be obtained. Moreover, if 
F and G are both linear in the scalar components of ^ and the Grassmann bracket i[F,G]c 
formed from scalar components of F and G will agree with the canonical anticommutator i{F, G}c 
formed from the corresponding operator scalar components, and will be a c-number. Thus, for 
linear F and G we can define a “Dirac anticommutator” {F,G}d by 


{F, G]d ={F, G}c - E E{^’ ^a}cT-;{^b, G}c 

a b 

Tab={^a,^b}c ■ 


(29) 


When one or both of F and G is bilinear, the Grassmann bracket i[F, G]c formed from the scalar 
components of F and G will agree with the canonical commutator formed from the corresponding 
operator scalar components, and we can define a “Dirac commutator” by a formula analogous 
to Eq. (I29h in which each anticommutator with at least one bilinear argument is replaced by a 
commutator. In this way we get a mapping of classical brackets into quantum anticommutators and 
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commutators, that inherits the algebraic properties of the Dirac bracket, including the chain rule, 
with the Jacobi identities for odd and even Grassmann variables mapping to the corresponding 
anticommutator and commutator Jacobi identities. 

To complete this correspondence, we must show that the Dirac anticommutator of and 
(with a = 1,2, /3 = 1,2 the spin indices, u = v = l,...,n the internal symmetry 

indices, and i = 1,2, 3, j = 1,2,3 the spatial vector indices) has the expected positivity properties 
of an operator anticommutator, by showing that for an arbitrary set of complex functions Af^{x), 
we have 


I dVyAr{x)Af\y){^ri^),^V"m}D > 0 . (30) 

We demonstrate this in several steps, in covariant radiation gauge. First we examine the conditions 
for positivity of the canonical anticommutator and Poisson bracket. 


J d^xd^yArix)Af^m^r{x),^]f^^mc = j d^xd^yAr{x)Af'^m^r{x),^y^mc 


(31) 


From = iPj^ — CjkiP^'^cr^^, we hnd that 


i^rix), ^Y\y)]c = - + iej,k^f)5-'’5Hx-y) 


1 


= - i{ajair^6^^6^{x - y) = -2i{5ij - -a^a^r^5^^5\x - y) 


(32) 


Multiplying by i/2, and writing + ilf^, i = 1, 2,3, a = 1, 2, tt = 1,..., n, with R and I 

real, the right hand side of Eq. m evaluates to (we suppress the internal symmetry index u from 
here on, so (Rf)'^ means > ®tc.) 

2=1 Q; = l 

If all three components Af, i = 1,...,3 are present, the expression in Eq. (1331) is not positive 
semidefinite. But when only two of the three components are present, as a result of application of 
a constraint, then each of the four squared terms on the right hand side of Eq. (|33p contains only 
two terms, and so the expression in Eq. (|33p is positive semidefinite by virtue of the inequality 


+ y2 _ 1 ± y)2 = i(x =F Yf > 0 


(34) 


IS 


Another way of seeing this, noted by both Velo and Zwanziger and Allcock and Hall 
that because (TjUj = 3, the expression Wij = 5ij — is not a projector. But when one 
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component of a, say 1 T 3 , is replaced by 0 , so that one has — Yli=i — ‘^1 then 


1 1 ^ ^ 

WiiWij = 6ij - 2-aiaj + -ai ^ aiaiaj 


— ^ij ^ O'id'i — 'Wt 




(35) 

i 1=1 ^ 

and Wij is a projector and hence is positive semidefinite. So we anticipate that proving positivity 
will require projection of Eq. (|32p into a subspace obeying at least one constraint on 

The next step is to use the property that the Dirac bracket of linear quantities F and G reduces 
to the canonical bracket of their projections into the subspace obeying the constraints, when (as 
is the case here) all constraints are second class, that is they all appear in the Dirac bracket 
Referring to Eq. (1281) . let us define 


a b 

G=G-EE[c.'f.icr„y-i-i, . 


a b 


(36) 


so that 


[F, Ck.lc =[F, Ck.lc - ^ ^[F, <^>a]cT-,\^b, ‘^clc 

a b 

= [F, ck.lc - ^ ^[F, ^a]cT-,^Tb, 

a b 

= [F,^c]c-J2^F,^a]cSac = 0 , 

a 

(37) 

and similarly for G. As a result of this relation, which holds when the canonical brackets are 
simply numbers (as in the case here where <hc and F, G are linear), together with symmetry of the 
canonical bracket [G, <kc]c = we see that 

[F,G],, = [F,G]c . (38) 

These properties of Eqs. (I36l) ~ (l38p carry over when we replace Grassmann numbers with opera¬ 
tors, and classical brackets with anticommutators, since in the linear case all anticommutators of 
linear quantities are c-numbers that commute with the operators, and since the anticommutator 
is symmetric. Thus we have 

mix),^y{y}D = mix),me 


( 39 ) 
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To further study the properties of Tj(x) and (with spinor indices suppressed), let us now 

return to our original labeling of the constraints by 4>a and Xa as in Eq. m, so that we have in 
the Dirac bracket formalism 

Ti(f) = Ti(x) - , (40) 

a b 

and a similar equation (with the roles of and Xa interchanged) for Tj(y), with a, b summed from 
3 to 4. We now note two important properties of this equation. The first is that it is invariant 
under replacement of the constraints Xa by any linear combination x'a — Xb^ba, with the matrix K 
nonsingular, since the factors K and K~^ cancel between x'a ^^4 . (More generally, the Dirac 

bracket is invariant under replacement of the constraints by any nonsingular linear combination 
of the constraints, reflecting the fact that the Dirac bracket is a projector onto the subspace 
obeying the constraints, and this subspace is invariant under replacement of the constraints by 
any nonsingular linear combination of the constraints.) The second is that if we act on Tj(x) with 
either or {a x D^)i, we get zero. For example, recalling that in covariant radiation gauge 
D^i'^i{x) = 04 (x), we have (with spatial variable labels x suppressed) 

Di^i = 04 - ^ Xa]cM~^ct)b = 04 - ^ ^ MiaM~^^(l)b = 04 “ ^ '54b0fe = 0 , (41) 

a b a b b 

and similarly for (cf x Dg)i, with 04 replaced by 03 . 

Let us now write 4'i(x) as a projector Rij{x,y) acting on T^^y), giving after an integration by 
parts on y, 

^iix) = j d^yRijix,y)^j{y) , 

Rij{x,y) =6ij5^{x - y) [ d? z[^i{x),Xa{z)]cM~^{z,y)'^bj{y) , 

a 6 4 

(42) 

with 

V3i(y) = (<? X , V4i(y) = ^yj ■ (43) 

By virtue of Eq. (I4ip and its analog for a x D, we have 

d^xi^ijix,y) 0 , 

{d X D^)iRij{x,y) =0 . 

( 44 ) 
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Since 


a ■ D^GiRij (f, y) = D^iRij(x, y) + i(a x Ds)iRij (x, y) , 
then assuming that a ■ D \s invertible Eqs. (1441) also imply that 


( 45 ) 


aiRij {x, y) = 0 . (46) 

Next let us focus on the bracket ['l'j(x), Xa(-?)]c appearing as the first factor inside the sum. Setting 
L = z3 in Eq. (fT8]l we have 

[^{x),X^{z)]c =‘^Ds 6 '^{x - z) , 

[^(x),X4(^)]c' ={iDx - (? X D^)S^{x-z) . 

(47) 

Using the invariance of 4'j, or equivalently of Rij, under replacement of xs, X 4 by any nondegenerate 
linear combination of xs? X4) l^t us choose the new combinations so that 

[^{x),X 3 {z)]c =(<? X D^) 6 ^{x- z) = ff 3 {x) 6 ^{x - z) , 

[^(x),X4(^)]c’ =Ds;S^{x-y) = f] 4 {x)d^{x - z) . 

(48) 


Substituting this into Eq. (j42p . we get the symmetric expression 

Rij{x,y) = 5ij5^{x-y) + '^'^ i d^zf]ai{x)M~^^{x,y)^bj{y) . (49) 

a b 

By virtue of this symmetry, the projector Rij is annihilated by the constraints D^j and {a x ^y)j 
acting from the right, which in turn implies that in addition to Eq. (I46p we also have 


Rij{x,y)aj = 0 . (50) 

An explicit construction of Rij{x, y) and verification of Eqs. ()46]l and ([501) is given in Appendix C. 

Returning now to Eqs. ()30p and ()39p . writing 'kf and in terms of projectors acting on 'kf 
and we have (using and continuing to suppress internal symmetry indices u,v, 
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which are contracted in the same pattern as the spatial vector and spin indices) 

I d^x I d^yAf{x)A*/{y){'^f{x),^]^mD 


= jd^xj d^yAf{x)A*/{y){^<t{x),^^{y)}c 

= J J j J d^w RfP{x, {z),^l^{w)}cR*^{y,w) 

^ I I j d^z j d^wRf[^{x,z)2(^6imd'^^5^{z - w)R*^{y,w) 

=2 j d^z J d^xAf{x)R^,^{x,^ J d^yA^{y)R'^^{y, ^ 


(51) 


which is positive semidefinite. 

We conclude that the anticommutator of ^ with is manifestly positive semidefinite in co¬ 
variant radiation gauge. The duality of the ^ 3^4 and X 3^4 constraints in this gauge is essential 
to reaching this conclusion; if gauge fixing were omitted, or if another gauge were chosen, this 
symmetry would not be present and we could not deduce positivity in a similar fashion. 


IV. ALTERNATIVE LAGRANGIAN AND HAMILTONIAN FOR THE T EQUATION 

IN COVARIANT RADIATION GAUGE 


Up to this point we have worked with the original action of Eq. (1-23) and the canonical 
momentum derived from it. We give here another approach, based on setting up an action for the 
T equation of motion from which Tq has been eliminated by the secondary constraint, 

HqT =DR • T zU X T , 

R={d- + d X E) , (52) 


which holds when the primary constraint x = 0 is obeyed. Consider the self-adjoint action 


S = j d^xL = ^ J d'^x^^ ■ {Do^-iD - DR-^ - R^D ■^) . (53) 

Varying with respect to and imposing two constraints: (i) the primary constraint x = d-D x T = 
0, and (ii) the gauge fixing constraint U • T = 0, we get the equation of motion of Eq. ([52]) . For 
the canonical momentum conjugate to T, we find 


P = 


d^S 

d{dQ^) 



(54) 
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which implies that 

= 2iP . (55) 


For the Hamiltonian corresponding to the new action, we find (again for simplicity taking Aq = 0, 
and integrating the middle term by parts) 


/' 


H = I Pxdn^ -P-L 


j cPx^^ ■ { - D X ^ - i^R ■ ^ + iR^3 ■ 

- ^ J Px^^ • z5 X ^ , 


(56) 


where in going from the second to the final line we have used the constraint 11 • ^ = 0 and its 
adjoint ^ = 0. The Hamiltonian H is again the energy integral calculated from the left chiral 
part of the stress-energy tensor, and expressed in terms of the canonical momentum is 

H = -i J PxP -Dx^ . (57) 

From here on the argument parallels that of Secs. 2 and 3, but is simpler. For the canonical 
bracket of Tj(x) with 'l'j(y) we have 

[Ti(x),Tj(y)]c = [^i{x),2iPj{y)]c = -2iSij6^{x - y) , (58) 

and so multiplying by i to convert to a canonical anticommutator we get 

{^i{x),^]{y)}c = 26ij6^{x - y) , (59) 

which is positive semidefinite. The complete set of constraints is 


4>3=X = ■ 3 x^ , 

(j)^ =3 ■ T , 

X3 =x^ = X ll • (T = 2iP • fj X ll, 

^ = 2iP • . 

(60) 


The constraints 4 >3 , are identical to ipi , ())2 of Eq. while the constraints X 3 ) X4 are xi, X 2 

of Eq. ()A1I) up to an invertible linear transformation (just interchange of the x constraints and 
division by 2 i). Thus the projector Rij{x, y) is the same as that calculated in Appendix A, and the 
Dirac anticommutator given by Eq. (I39h is is positive semidefinite by Eq. (I51h . this time without 
using the fact that Rij is projected to zero by ai and aj. 
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V. LORENTZ COVARIANCE OF COVARIANT RADIATION GAUGE AND LORENTZ 

INVARIANGE OF THE DIRAC BRACKET 

We study next the behavior of covariant radiation gauge and the Dirac bracket under Lorentz 
boosts. The Rarita-Schwinger field and its left-handed chiral projection both have a four- 
vector index ^ and a spinor index a. Under an infinitesimal Lorentz transformation, the trans¬ 
formations acting on these two types of indices are additive, and so can be considered separately. 
The spinor indices are transformed as in the usual spin ^ Dirac equation by a matrix constructed 
from the Dirac gamma matrices, which commutes with Hence the spinor index transformation 
leaves the covariant radiation gauge condition D ■ T invariant. 

This leaves the transformation on the vector index to be considered, and this is a direct analog of 
the Lorentz transformation of radiation gauge in quantum electrodynamics . Since the radiation 
gauge condition is invariant under spatial rotations, we only have to consider a Lorentz boost, 

X -^x' = X + vt , 
x^ =t ^ t' = t + V ■ X . 

(61) 

Under this boost, the field T transforms as 

T T' = T -b . (62) 

For an observer in the boosted frame, covariant radiation gauge would be D^/ • T' = 0, with 
= VS' + where A' = A + 0{v). Applying this to T'(T',t') and using the covariant 
radiation gauge condition in the initial frame, we get 

Ds' ■ T' = Vj^j{x,t) , (63) 

with 'Sj{x,t) a local polynomial in 'L, Tq and the gauge fields, where we have dropped primes on 
the right hand side since there is an explicit factor of v. So in the boosted frame T' does not obey 
the covariant radiation gauge condition, but this can be restored by making a gauge transformation 

T' ^ - D{D^)-^VjJ:j{x, t) . (64) 

Hence the covariant radiation gauge condition is Lorentz boost covariant, although not Lorentz 
boost invariant. 
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Referring now to Eq. (lAlOh . we note that the covariant radiation gauge Dirac bracket and 
the anticommutation relations are invariant under infinitesimal Rarita-Schwinger gauge transfor¬ 
mations, such as that of Eq. (f64|l . up to a remainder that is quadratic in the gauge parameter. 
Hence the covariant radiation gauge Dirac bracket and the anticommutation relations following 
from it are Lorentz invariant, since a finite Lorentz transformation can be built up from a series of 
infinitesimal ones. 


VI. PATH INTEGRAL QUANTIZATION 


An alternative method of quantization to the Dirac bracket approach is setting up a Feynman 
path integral. Again, we will specialize to the case where the external gauge potentials, and hence 
D, are time independent, since the simplest discussions of path integrals for constrained systems 
assume time-independent constraints. As noted above, this assumption can be dropped when the 
gauge field is quantized along with the Rarita-Schwinger field, leading to a more complex system 
of constraints and constraint brackets. 

When the constraints are time independentjthe classical brackets of Eqs. (1-57) and (1-58) have 


the form needed to apply the Faddeev-Popov method for path integral quantization. (This has 
been applied in the free Rarita-Schwinger case by Das and Freedman^] and by Senjanovic 


The general formula of 
factor) reads 


m-) 

8|] for the in to out S matrix element (up to a constant proportionality 


(out|5|in) oc f exp (^iS{q,p)) n dp{q{t),p{t)) , 

d t 

dp{q,p) = n SiXa)d{(j)a)idet[(f>a,Xb])^ n dpidqi , 


(65) 


where ^ = 1 when all canonical variables are bosonic, and ^ = — 1 in our case in which all canonical 
variables are fermionic, or Grassmann odd. In applying this formula, we note that since the action 
S of Eq. (T23) and the bracket matrix M of Eqs. (I-59)-(I-62) are independent of and P^u 
we can immediately integrate out the delta functions in these two constraints. Also, since the 
canonical momentum P is related to by the constant numerical transformation of Eq. (T54), 
we can take as the integration variable instead of P, up to an overall proportionality constant. 
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So we have the formula, after an integration by parts in the second term, 

(out|5'|in) oc J exp (i- J '1'qit • D x ^ ^ ^ x + — 'h^-ux Hq^]) 


t.X 


= J exp J d^rc[—'P qX — x^'I'o + — ^’^•ax Hq^]) 

x]Jo!;u(«'o,^cu^,^'^) • 

t,X 


( 66 ) 

Here 

= ^JJ<J(Xc)5(<('c)^ (det[())a,Xb])"^(i^orf^o'^^'^^^ , (67) 

with d'ho and d'l'Q each a product over the spinor components, and and each a product 
over the spinor-vector components. 

As our next step, we can carry out the integrations over 'hg and 'hj, using the delta functions 
5 {4>2) and S{x 2 ) ■ This leaves the formula 

(outlAlin) oc J exp (i^ J d^x[—• (H -|- cj x .B)(it • .B)“^X 

-xH^ ■ + dxE)-'^ + '^^-Dx^-'^^- ax HqT]) 

xnd^(T,#) , 

t,X 

( 68 ) 


with 

d^(T,T'f) = ^]J(5(xc)d(0c)^ (det[(/>a,Xfe])"^d^d^^ , (69) 

so that only the remaining constraints </> 3 , 4 , X3,4 are used in constructing the determinant 
det[(l)a,Xb]- 

Finally, using the delta functions 5 {(j)3) = (5(x) and ^(xs) = d(x^) to simplify the exponent, we 
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end up with the elegant formula 

(out|5|in) oc J exp J • [Zl x ^ x Hq^]) 


t,X 


(70) 


which as in Dirac bracket quantization, employs as Hamiltonian the energy integral computed in 


Eq. (1-35) from the stress-energy tensor. In using this formula, the customary procedure 


111 ] would 


be to put the bracket matrix that is the argument of the determinant back into the exponent by 
introducing bosonic ghost fields (pc- 

VII. CONCLUSION AND DISCUSSION 


To conclude, we see that when a covariant radiation gauge constraint is included, the problems 

n n 

with canonical quantization found in and [J| are avoided: The Dirac bracket is well-defined in 
the limit of zero external fields, and is positive semidefinite. Thus our conclusion in (I) that the 
classical theory of gauged Rarita-Schwinger fields is consistent extends to the quantized theory of 
gauged Rarita-Schwinger fields as well. As noted in (I), this means that in constructing grand 
unihed theories, one can contemplate an anomaly cancellation mechanism in which the gauge 
anomalies of Rarita-Schwinger fields cancel against those of spin-^ fields, as first suggested in 12l | 
and as used in the SU{8) family unification model of {l^ . 

Some final remarks: 

1. In quantizing, we assumed that the gauge fields A are time independent, so that d/dt and 
D commute. As noted, this assumption can be dropped if the gauge fields are treated as 
dynamical variables, leading to an extension of the bracket structure, involving fermionic 
brackets as well as bosonic ones. (For a discussion of bosonic versus fermionic constraints, 
see 


Q. 


2. In demonstrating positivity of the anticommutator in Sec. 3 (but not in Sec. 4), we used the 
condition u • T = 0. Deriving this from the covariant radiation gauge condition D • T = 0 
assumed the invertibility oi a ■ D, and attainability of covariant radiation gauge assumed 
the invertibility of (D)^. The conditions for invertib ility of these two o per ators remain to be 
studied. (The open space index theorems of Callias 15| and Weinberg 15| involve a ■ D + i(l>, 


with (p a scalar field, and so do not give information about the invertibility of a ■ D.) 
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Appendix A: Construction of the projector Rij{x,y) 


Since there are only two (pa constraints and two Xa constraints, we index them a = 1,2 rather 
than o = 3,4 as in the text, and use the invariance of Rij{x,y) under changing the linear combi¬ 
nation of the Xa constraints. We start from the constraint set 

(pi = a X D , xi = P , 

<p2 = D - ^ , X2 = P ■ ■ 

(Al) 


For the bracket matrix 


Mab{x,y) = [4>a{x),Xb{y)]c 



we find the matrix elements 


(A2) 


A = — iga ■ B5^{x — y) , 

P ={^{Dxf + 9^ ■ B)S^{x -y) = S^{x - y){2{5gf + ga ■ B) 
C ={Dx?S^ix -y) = SHx - , 

V =iga ■ BS^{x — y) . 


(A3) 
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We write the inverse matrix M ^{z,w) as 



which obeys 



(A4) 


(A5) 


In terms of the inverse matrix, the projector Rij{x,w) is given by (with internal symmetry indices 
suppressed) 


Rij{x, w) =6ij6^{x — w)l 

'?c)(cJ X DxiQi^X yj j 

+{a X D^)i'H{x -w){ax + (ct x D^)iZ{x - w)^^j 


(A6) 


From this expression, we find 

D^iRij{x,w) = Rij{x,w)^^j = {a X D^)iRij{x,w) = Rij{x,w){a x = 0 . (A7) 

In verifying these, it is not necessary to evaluate the inverse matrix; instead, after contracting on 
the vector index i or j one expresses the resulting pre- or post- factor in terms of and 

then uses the algebraic relations following from multiplying out the matrices in Eq. (|A5p . Finally, 
contracting 

(7 • (yZ)x A ^ ^x')i ; 

* ^5^ (^5^ i(T X , 

(A8) 


with Rij{x,w), we conclude that 

aiRij{x,y) = Rij{x,y)crj = 0 , (A9) 

when ff • H is invertible. 

As a consequence of Eqs. ()42ll and ()A7p . 'l'j(x) is invariant under the transformations 

+ , 

^ + a X De . 


(AlO) 
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The first of these implies that the canonical anticommutation relations are invariant under in¬ 
finitesimal Rarita-Schwinger gauge transformations starting from covariant radiation gauge. 
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